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Exchange of dominance between twin kHz quasi-periodic oscillations (QPOs) observed in some low-mass-X-ray-binaries (LMXB) 
suggests the possibility of a resonance between two oscillatory modes. We study the behaviour of the effective gravitational potential 
around specific resonant radii, and estimate the role of the higher-order terms governing the non-linear, anharmonic forcing. We 
discuss the impact it has on the mode amplitude in the linear and non-linear regimes. We also discuss a related possibility of lowering 
of the neutron star mass estimates from the highest observed QPO frequencies. 
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1. Introduction 

High-frequency quasiperiodic oscillations (HF QPOs) of the X- 
ray brightness have been observed in several accreting binary 
systems containing neutron stars (NSs, see van der Klis 2000, 
2006; Belloni et al. 2005) and black holes (BHs, see McClintock 
& Remillard 2004; Remillard & McClintock 2006). In NSs they 
often appear in the form of two distinct modes (so-called "twin 
peaks") with frequencies correlated with the X-ray intensity and 
covering a relatively wide range ~ 50 - 1300 Hz. Henceforth, 
we adopt the convention of referring to twin-peak QPO modes 
as to lower and upper QPOs and denote their frequencies as vl 
and vrj. 

The ratio of HF QPO frequencies when observed in black- 
hole systems in pairs is usually exactly, or almost exactly 3 : 2, 
while in the NS systems the ratio R = v v lv L is concentrated 
around R — 3 : 2 (e.g., Torok 2005; Abramowicz et al. 2005). 
Most HF QPO models involve orbital motion in the inner re- 
gions of an accretion disc. Because of the 3 : 2 ratio appearance, 
several proposed models consider resonances between the QPO 
modes (Abramowicz & Kluzniak 2001 ; Abramowicz et al. 2003; 
Lamb & Coleman 2003; Kluzniak et al. 2004; Petri 2005a,b; 
Torok et al. 2005; Sramkova et al. 2007; Stuchlfk et al. 2007, 
2008a; Stuchlfk & Kotrlova 2009; Mukhopadhyay 2009, and 
others). 

Quite recently Torok (2009) discovered that, in the six atoll 
NS systems, the rms amplitude difference A r = r^-r\j is zero for 
the resonant frequency ratios R = 3 : 2 at a frequency v^ 2 specific 
to a particular system. The quantity A r is positive for frequencies 
above v^ :2 and negative for frequencies below v^ 2 . 1 Similar "en- 
ergy switch effect" also occurs at the frequencies corresponding 
to the resonant ratios 4:3 or 5 : 4. Nevertheless, these ratios cor- 



1 It has been argued by Horak et al. (2009) that a natural explanation 
of such an "energy switch" effect arises in the framework of the non- 
linear resonant phenomena. 



respond to very high frequencies where kHz QPOs disappear. 
The "disappearance" occurs since the QPO amplitudes and co- 
herence times strongly decrease; see Barret et al. (2005a), Barret 
et al. (2005b, 2006), Mendez (2006), Torok (2009), Torok et al. 
(2008a,c,b), and Boutelier et al. (2010). Figure 1 illustrates both 
the energy switch effect and the QPO disappearance. 



We suggest that the inner edge of the disc can be related to 
the resonant radius corresponding to the 5 : 4 (4 : 3) frequency ra- 
tio. We study properties of the effective potential of geodesic mo- 
tion in the vicinity of the resonant points corresponding to these 
frequency ratios. For identification of the QPO frequencies we 
assume a widely discussed relativistic precession model (hence- 
forth RP model) where vl = vk - v r and vu = vk with vk and 
v, the Keplerian and radial epicyclic frequency of the geodesic 
motion (Stella & Vietri 1998, 1999; Karas 1999; Kotrlova et al. 
2008). This identification roughly merges with a model assum- 
ing m=-l radial and m=-2 vertical disc-oscillation modes (see 
discussion in Torok et al. 2007, 2010). We expect that our anal- 
ysis is relevant not only to these two specific models but also 
to some others that assume resonances between the two QPO 
modes in the regions close to the innermost stable orbit and a 
frequency ratio decreasing with decreasing orbital radius. 



We focus attention on the potential depth at the resonant radii 
and related maximal amplitude of the oscillations, and we deter- 
mine the energy levels for the harmonic oscillations given by the 
expansion of the potential around the resonant radii (whereas 
the second order term of the expansion gives the frequency of 
the epicyclic motion). Furthermore, we estimate the role of the 
higher order expansion terms of the potential that govern the 
non-linear terms relevant to the resonant phenomena (Nayfeh & 
Mook 1979). 
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Fig. 1. Energy switch effect (black datapoints and lines from Torok 2009, see this paper for details) and the twin QPOs fadeaway at low-frequency 
ratios. The colour-lines after Barret et al. (2005a,b, 2006), Mendez (2006), and Torok (2009) indicate the drop in the behaviour of the total 
energy of the two QPOs (red and green lines) and the drop in the high cohererence of the lower QPO (orange and yellow lines), a) Frequency 
ratios emphasized in 4U 1636-53 and 4U 1608-52 (3:2 and 5:4) respectively 4U 0614+09 and 4U 1728-34 (3:2). b) 3:2 and 4:3 frequency ratios 
emphasized in the two sources 4U 1820-30 and 4U 1735-44. 



2. Effective potential near the resonant radii at the 
Schwarzschild spacetime 

Using the dimensionless Schwarzschild radial coordinate x = 
r/M, the effective potential of the radial motion in the 
Schwarzschild spacetimes can be given in the form (Misner et al. 
1973) 



V*(x;L) = 1-- 1 + , L\x) 



x-3 



(1) 



where L is the specific angular momentum of a test particle eval- 
uated for the specific case of a circular orbit of given r (x). 

In the Schwarzschild spacetime, the orbital Keplerian fre- 
quency and the radial epicyclic frequency read (Aliev & Galtsov 
1981; Torok & Stuchlflc 2005) as 



F 6 

V K = -T7T, V r = V K 1 

X ilL \ X 



1/2 



F = 



2nGM 

The resonant radii for the RP model are determined by 

6n 2 



' 7 K 



Vk — Vr "l 



m (2n — m) 



(2) 



(3) 



We thus find x-}a = 6.75, X4.3 = 6.40, andx^ = 6.25. The 
relevant effective potentials are then given by 



y eff (x; L) = V eS (x; L = L n:m ), 



L . = 



These potentials can be expressed in the form 

/ 2\l I2n 4 x- 2 
V eS (x;L n:m ) = 1 - - 1 + — — — -5- 

\ xj\ m(2n — m)(ln^ — Inm + m z ) 



(4) 



(5) 



and are presented for n : m = 3 : 2, 4 : 3, 5 : 4 in Fig. 2. Then we 
can determine the potential depth at the resonant radii giving the 
energy barrier for the stability of the circular motion, 



^F n - m — V ma x(-*-5 Ln:m) ^min(-*-? ^n:ni}- 



(6) 



The results are given in Table 1. Similarly, we can determine 
the maximal amplitude A n:m allowed for the oscillations in the 
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Fig. 2. The effective potentials L n:m determined from the RP model 
(black lines). The blue dashed lines denote the Taylor approximation 
up to the quadratic terms. The red solid lines correspond to the case 
of A p /A r ~ 0.97 and the red dashed lines correspond to the case of 
A p /A r = 0.99. 



vicinity of the resonant radii (see again Table 1). We also give the 
ratio - 221 , which demonstrates the maximal relative magnitude of 
the oscillations at the resonant radii. We can see that the potential 
depth and the maximal amplitude at the resonant radii strongly 
decrease with frequency ratio decreasing from 3 : 2 down to 5 : 4. 
Nevertheless, in all three cases, the maximal amplitude of the 
oscillations remains relatively large. The relative amplitude is 
largest for the 3 : 2 ratio when ~ 0.34. We note that the 
potential depth at ^5 4 is by more than one order (~ 1 /20) smaller 
than those at X3 : 2, while the maximal amplitude of the oscillatory 
motion is smaller only by factor of ~ 3. 

The a; coefficients in Taylor expansions of the effective 
potential about the resonant radii are 
for X3-2'- 

3.4444xl0~ 4 , a* = 3.40148 xlO" 5 , 



ai =0.944118, 
-3.7857xl0~ 5 , 
8.32113xl0~ 7 , 
for 

ai = 0.943242, 

-6.97262xl0~ 5 

1.53105xl0" 6 , 



a 2 = J.4444X1U , a 3 : 
a 5 = 1. 34255 xl0~ 5 , a 6 : 
a 8 = -1. 78353 xlO" 7 ,.. 



-3.58131xl0~ 6 , 



a 2 =2.37896xl0- 4 , 
a 5 =2.35669xl0- 5 , 
a 8 = -3.43018X10- 7 



a 3 = 1.11514xl0- 4 , 
a 6 = -6.36989xl0- 6 , 



a 4 : 

;i7 : 



a 4 = 
a 7 = 
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Table 1. The potential depth AE nm at the resonant radii given by 
relation (6), the maximal amplitude A mm and the relative amplitude 
A nm jx nm allowed for the oscillations in the vicinity of the resonant radii 
for the RP model. 



mm 


3:2 


4:3 


5:4 


AF 


3.27 x irr 4 


5.77 x irr 5 


1.51 x 10- 5 


A 


2.31429 


1.21008 


0.75251 




0.34286 


0.18908 


0.12040 



Table 2. The potential depth AE( Vp -v MlN ) and the related maximal am- 
plitude A p of the harmonic oscillations for two precision levels given by 
the ratios A p /A,. = 0.99 and A p /A r = 0.97. 



n : m 


3:2 


4:3 


5:4 


A p /A r > 0.99 


AF 

A £(fp-fMIN> 


9.5 x 10~ 6 
34.4 

0.166085 


3.9 x 10~ 7 

148.1 
0.0404892 


6 x 10" 8 
251.8 

0.0189513 


A, 


0.167759 


0.0408933 


0.0191297 


A p /A r > 0.97 


A£ (v p -v M in) 


45 x 10" 6 


8 x 10" 6 


3 x 10" 6 




0.722944 


0.366759 


0.268010 


A, 


0.729074 


0.370491 


0.272475 



for X5 : 4! 

ai = 0.94299, a 2 = 
-8.98277xl0~ 5 , a 5 
2.00698 xl(T 6 , a 8 = 



1.67063 xl0~ 4 , 
= 3.00803xl0- 5 , 
= -4.59198x10" 



a 3 : 



1. 60380 xl0~ 4 , 
-8.20919xl0~ 6 , 



a 4 : 
a 7 : 



The expansion coefficients clearly demonstrate the importance 
of the higher order terms. For the sake of completeness we give 
the expansion terms up to the eigth order, while for the non- 
linear resonant phenomena usually the third and fourth terms 
are relevant (Nayfeh & Mook 1979; Stuchlfk et al. 2008b). We 
can see that at the x\2 radius the harmonic (second order) term 
is by one order higher in comparison with the third and fourth 
order terms. On the other hand, at the x^ radius all the second, 
third, and fourth order terms are of nearly the same magnitude, 
indicating a strong role of the non-linear resonant phenomena 
even for small perturbations from the purely circular motion. At 
the JC4 : 3 radius, all three terms are of comparable magnitude. 

Now we can determine the regions that extend around the 
resonant radii where the oscillation can have a harmonic charac- 
ter, i.e., where only the second-order term of the Taylor expan- 
sion of the effective potential is relevant. We also determine the 
maximal energy level and the maximal amplitude of harmonic 
oscillations. We carry out the estimate for two precision levels 
given in terms of the ratio of the amplitude governed by the total 
effective potential (A r ) and the amplitude implied by the sec- 
ond order (harmonic) part of its expansion (A p ). The results are 
presented in Table 2. Clearly, both the energy level of the har- 
monic oscillations and their amplitude strongly decrease with 
the resonant radius (and frequency ratio) decreasing. Table 2 
gives the evidence for an increasing role of the non-linear res- 
onant phenomena when the resonant radius decreases. The har- 
monic character of the oscillations near x^\ is only possible for 
AE(v p -v Mm ) ~ 6 x 10~ 8 , which is by factor ~ 250 smaller than 
the related potential depth. 



RP model 3:2 
TO model 4:3 
5:4 
1:1 




Fig. 3. The x„ „, (j; n, m) functions for the RP model (red lines) and the 
TO model (see Sect. 4). The thick green line represents the radius of the 
marginally stable orbit x ms where v v /v L = 1 for both models. 

3. Resonant radii at the Kerr geometry 

We next assume Kerr geometry describing the gravitational field 
around rotating BHs (Bardeen et al. 1972). This is in general a 
somewhat worse, but applicable, elegant, and fully analytical ap- 
proximation to the case of the more complicated Hartle-Thorne 
geometry, which describes rotating NSs quite well (Hartle & 
Thorne 1968; Berti et al. 2005). Moreover, for high NS masses, 
such as those implied by the RP model, the Kerr geometry ap- 
proximation is very good (Torok et al. 2010). 

Using the dimensionless Boyer-Lindquist radial coordinate 
x = r/M and the dimensionless spin j the effective potential of 
the radial motion in the Kerr spacetimes can be given in the form 



V eS (x; j; L) = 



(7) 



2jL + V j 2 -2x + x 2 y/x (2 j 2 + fx + xV- + x 3 ) 
x 3 + j 2 (2 + x) ' 

where L denotes the specific angular momentum of a test par- 
ticle. When the particle orbits in a corotating circular orbit at 
given r (x), its specific angular momentum reads as 



-2jx 1 ' 2 +j 



Ltx;j) = 



1/2 



(8) 



2>x + 2 jx l/1 

In the Kerr spacetime the orbital Keplerian frequency and the 
radial epicyclic frequency read as 

(9) 



v K = 



+ ./ 



6 8; 



3/ 



Again, the zero point of the radial epicyclic frequency defines 
the marginal stable geodesies x ms that is implicitly given by 

Ax 112 1 / 

j = Us(x) = — ~ 3 V3* 2 - 2x. (10) 

The resonant radii for the RP model are now implicitly given by 



J = Jn:m U; n,m) = — 



4- J3x 



1 



(n — m) 2 



(11) 



and presented in the explicit form x n:m (j) in Fig. 3. The angular 
momentum of test particles orbiting at the resonant radii is given 
by 



2jx l J 2 (j) + j 2 



x n:mU) ^ Xn:m(j) ' 



U) + 2jx x J. 2 {j) 



(12) 
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Fig. 4. Effective potentials for the RP and TO models. The green curve L ms corresponds to R = 1 for both models. Also RP3 : 2 = TO5 4. 
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Fig. 5. a) The spin j profile for the potential depth giving the energy barrier of the stable circular motion at the resonant radii AE„-„,(j) for the RP 
model (top panel). The relative potential barriers AEi-^ij)/ AE„ ln (j) are illustrated for the considered resonant points (bottom panel), b) The same 
as a) but for high-spin geometries, c) The maximal amplitude A„ m (j) of oscillations about the resonant radii (top panel) and the ratio -^(j) for 
the RP model (bottom panel). 



The effective potential related to the angular momentum L n:m 
reads as 

V eS (x;j;L n . m (x n . m ; j)) = (13) 
2jL„ :m + V;' 2 - 2x + x 2 -^x {ifi + fix + xL; rm + x 3 ) 
x 3 + fi (2 + x) 

For selected values of the spin (j = 0, 0.4, 0.95) the effective 
potential with minimum at the resonant points is given in Fig. 4. 
In Fig. 5 we give the spin j profile for the potential depth at the 
resonant radii AE n:m (j), the related maximal amplitude A„ :m (j) 
of oscillations about the resonant radii, and the ratio -f^ij). 

We can see that the potential depth AE n:m (j) and the ratio 
-^(j) with growing j increase. This increase is most restricted 
for the resonant radius xs : 4 and is greatest for xj,a. The amplitude 
of the oscillations A„ :m decreases when j grows up to ~ 0.9, and 
then it increases. On the other hand, the relative magnitude of 
the oscillations grows within the whole allowed interval of 
j. This indicates that the potential detectability of the resonant 
oscillations grows with increasing j. Especially near X3 : 2 the rel- 
ative amplitude is very large, > 0.34. 



4. Discussion and conclusions 

The role of the non-linear terms determining the non-linear res- 
onant phenomena grows for the decreasing resonant radius cor- 
reponding for the RP model to a decreasing frequency ratio. 
Extension of the possible harmonic radial oscillations decreases 
with decreasing resonant radius. This decline is strongest at xs : 4, 
reaching almost two orders, which is much higher compared to 
those at X3 : 2 where it is higher by only one order. The expected 
magnitude of oscillations around the resonant radii can then be 
quite large when the non-linear regime is entered, and it can be 
of the same order as the resonant radius. We note that the observ- 
ability of such phenomena is expected even for oscillations with 
the amplitude by one order smaller than x„ M (Bursa et al. 2005; 
Schnittman & Rezzolla 2006; Abramowicz et al. 2007; Bursa 
2008). 

If the edge of the accretion disc is located at the resonant 
radius x n:m (with n : m = 4 : 3 or 5 : 4) rather than at the 
marginally stable orbit x ms , the maximal frequency v max ( ) ob- 
served in a given source has to be attributed to the Keplerian 
frequency vk(x„ :) „) rather than to VK(x ms ), which is often applied 
in QPO-ISCO-estimates for the mass of the NS (e.g. Barret et al. 
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Table 3. Quantities A£„ :m , A nm , and A n:m /x n:m for the TO model. 



Fig. 6. Lowering of the estimated neutron-star mass. 



2005a,b, 2006). Therefore, we expect lowering of the NS mass 
estimated from the highest QPO frequencies owing to the shift to 
higher radii of the orbital motion. Since v max = v^(x ms (M ms )) = 
VK(xn:m(M n : m )), we immediately arrive at the relation 



M„ 



, 3 / 2 



+ J 



M„ 



v 3/2 



+ J 



(14) 



which determines the lowering of the estimated NS mass (see 
Fig. 6). Assuming the Schwarzschild geometry (j = 0), we find 
the relations 



M- 



5:4 



0.9406 M„ 



M 4:3 = 0.9077 M„ 



(15) 



The growing spin causes a further shift, and at the edge of the 
spin interval for the slowly rotating NS we find 

MsaU ~ O- 4 ) = 0.9340 M ms , M 4 - 3 (j ~ 0.4) = 0.8977M ms . (16) 

Such corrections are relevant and tend to lower the estimated 
masses down to the canonical value M ~ 1 .4 M Q . 

Throughout the paper we have assumed the RP model fre- 
quency relations. For comparison we discuss the behaviour of 
the effective potential in the Schwarzschild spacetime at the res- 
onant radii relevant for an alternative kHz QPO model, namely 
the model of deformed disc oscillations trapped in the gravita- 
tional field of the compact object that was introduced by Kato 
(2009) (henceforth the TO model). In the TO model the reso- 
nance condition reads (2vk - v,-)/2(vk - v r ) = n/m. We show 
the effective potentials with minima at the resonant radii of the 
TO model for the Schwarzschild spacetime, and for selected val- 
ues of spin j in Fig. 4. For the Schwarzschild spacetime we 
give the resonant radii, the potential depth, and the maximal am- 
plitude of the oscillations at these radii in Table 3. Notice that 
= x 5 4- Therefore, for all spin values, the resonant radii, po- 
tential depth and maximal amplitude corresponding to the TO 
model are greater than those related to the RP model. 
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